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We consider a nonlinear, strongly coupled, parabolic system arising in the modeling of burglary 
in residential areas. The system is of chemotaxis-type and involves a logarithmic sensivity function 
and specific interaction and relaxation terms. Under suitable assumptions on the data of the 
problem, we give a rigorous proof of the existence of a global and bounded, classical solution, 
thereby solving a problem left open in previous work on this model. Our proofs are based on the 
construction of approximate entropies and on the use of various functional inequalities. We also 
provide explicit numerical conditions for global existence when the domain in a square, including 
concrete cases involving values of the parameters which are expected to be physically relevant. 

1. Introduction 

1.1. Model and main results. In a series of recent papers, models based on partial differential 
equations have been derived to study crime, see [3] , [HI HH EU [22] . Some of these papers are 
most related to modeling burglary of houses. A basic issue here is to obtain patterns that describe 
the location of hotspots, and study phenomena such as appearance and disappearance of them, 
their stability and movement. 

A very successful model was obtained recently by Short et al. in 22 . They first derived an 
agent-based statistical model to study the dynamics of hotspots, taking two sociological effects into 
account: the 'broken window effect' and the 'repeat near-repeat effect'. The first effect refers to 
the observation that crime in an area leads to more crime, and the second one to the observation 
that houses burglarized at some moment have an increased probability of being burglarized again 
for some period of time after that moment. 

The agent-based model considered relies on the assumption that criminal agents are walking 
randomly on a two-dimensional lattice and are committing burglaries when encountering an attrac- 
tive opportunity. An attractiveness value is assigned to every house, which measures how easily 
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the house can be burgled without consequences for the criminal agent. In addition to be walking 
randomly, the criminal agents move toward areas of high attractiveness values. In turn, when 
a burglary occurs, it increases the attractiveness of the house that was burglarized and of those 
nearby. If no additional burglaries occur, then the local attractiveness decays toward a constant 
value. 

In a second step, by taking a suitable limit of the equations for the discrete model, the authors 
in [22) obtained the following system of parabolic differential equations (in adimensionalized form): 

( r)A 

— = nAA + NA + A -A. 
dt 

< (1) 
dN - f — 1 NA + A-A*. 



I dt V 



Here and thoughout this article, we denote 

0(A) = xlogA, with X > constant. (2) 

The model described in |22) involves the specific value x = 2 (we note that the parameter \ cannot 
be scaled out by a linear change of dependent or independent variables). The functions A(x;t) 
and N(x, t) respectively represent the attractiveness value and the criminal density at position x 
and time t. The first equation describes the evolution of the attractiveness of individual houses 
to burglary, and the second describes the burglar movement. Here 17 is a bounded domain in R 2 , 
r\ > is the diffusion rate of attractiveness, A is the intrinsic attractiveness, A is a contant that 
in the equilibrium case represents the average attractiveness. 

A first effort to study the dynamics of the model represented by problem ((TJ) was done in 
[19) . where a corresponding initial-boundary value problem with no flux boundary conditions is 
considered. To deal with this problem, the authors in 19 assume, for simplicity, that fl is a square 
and under some symmetry conditions are able to map this problem into one with periodic problem 
boundary conditions. Their main result is local existence of a solution, but global existence is 
left open. In addition, and as a simplification of the model in [22] . they considered a generalized 
version of a Keller-Segel chemotaxis model with the goal of understanding possible conditions for 
global existence vs blow-up of solutions in finite time for the original model. 

In [22], an important hypothesis is that burglars are generated in the model at a constant rate 
and leave the lattice immediately after they have committed a burglary. In [18J . a modified agent- 
based model for burglary was obtained, a new condition is introduced to modify these effects. This 
new condition, that is refer as burglar fatigue, models the effect that if burglars are sufficiently 
deterred, they will eventually get tired and will stop looking for houses to burglar. This consider- 
ation introduces changes to the original equations of Short et al. Some additional changes to the 
original equations come in by considering in the model that, if too many burglaries occur at some 
location, then burglars are likely to assume that most of the high- value is gone or that the owners 
of the place have implemented more strict security measures. 

As in the original model, by taking the continuous limit to this new agent-based model, a new 
system of parabolic equations is obtained in [18] . This system, complemented with initial data 
and no-flux boundary conditions, leads us to study the following initial-boundary value problem 
(again in adimensionalized form): 
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— = V AA + tpNA(l 



dN 
~dt 



A) 

VN - iVW(A) 



A- A, x £ ft, t > 0, 



dA dN 
ou ou 

A(x,0) = A (x), N(x,0) = N (x), 



X £ CI, t > 0, 



x£dfl, t>0, 



£ G ft. 



(3) 



Here again the function i? is defined by ([2]) with x = 2 in [15] . Moreover, it is assumed - as we will 
do throughout this paper - that 

ft is a sufficiently smooth bounded domain of R 2 , or a square, (4) 

rj, ip, uj, A > are constants (5) 

and that the initial data satisfy 

(A , N ) £ H 1+ P(ri) x L 2 (Q) for some p > 0, with A > in H and N > a.e. in ft. (6) 

The outward normal vector on 9ft is denoted by v. (We point out that all the functional analytic 
properties needed for the proofs remain true in a square, although the latter is not smooth at the 
corner points.) Problem ([3]) is locally well posed, the positivity of A being of course understood 
as part of the definition of solution (see Section 2 for details). 

In [Tg], some results about linearized stability /instability of the homogeneous steady-states 
and some numerical simulations suggesting the existence of hotspots were given for system ©. 
However, the (local and) global existence of solutions was left open. 

In this paper, our main goal is to give a rigorous proof of the global existence and boundedness 
of solutions of system ([3]) under suitable assumptions. In what follows, we denote 

£ = £o(ft) = H~ 2 K- 1/2 , 

where [i is the best constant in the Poincare-Sobolev inequality ([28]) and K is the best constant in 
the interpolation estimate ([50]) . Also, recalling that Aq £ C(ft) by Sobolev imbedding, we set 



max{l, A, max Aq} > A n 

n 



min{l, A, min Aq} > 0. 

n 



Our main result is the following: 
Theorem 1.1. Assume (0), (^-(0) and 

'A \2 



max -^-min )max{||iVo||i,|ft|} < £ r/^ \ 2 ■ 

Then the solution of problem {3|) is global and satisfies the uniform bounds 

sup ||^4(i)|| oo < oo 
t>o 

and 



sup || AT(t) Hoc < oo, for all r > 0. 



t>T 



(7) 

(8) 
(9) 



On the other hand, when \ < 1, by using an argument from [5], we can show that the conclu- 
sion of Theorem 11.11 remains true under a simple assumption on Aq and A and without any size 
restriction on Nq. However, this does not apply to the criminological model in [TJ] where x = 2. 
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Theorem 1.2. Assume |I|), (Rp-(0|), with < % < 1, A < 1 and ma,xAo < 1. Then the solution 
of problem is global and satisfies the uniform bounds 



1.2. Explicit global existence conditions and discussion. Interestingly, in view of practical 
applications, the constant appearing in Theorem 1 1 . 1 1 can be estimated explicitly in the case of a 
square. 



Theorem 1.3. Let fi = (0, L) , with L > 0. Then the result of Theorem \l.l \ is true with 

With Theorem 11.31 at hand, we shall now explicitly compute our sufficient condition for global 
existence and boundedness in concrete cases involving values of the parameters which are expected 
to be physically relevant (cf. |22[[T8] ). The parameters of system ([3]) are expressed in [18] in terms 
of parameters of the dimensional form of the model. To avoid confusion, dimensional parameters 
are marked with a hat in what follows. 

The domain is taken to be a square lattice where exactly one house is located at each lattice 
site. Distances are thus measured in units of house separations and the adimcnsional side length 
of the domain is thus essentially equal to ^/p, where p 3> 1 is the total number of houses. After 
a suitable renormalization (involving the diffusion and mean lifetime parameters of the burglars - 
see table 1 and formula (3.3) in |18[ p. 406]), the measure of the adimensionalized spatial domain 
Q. is given by \ft\ = p/3500. We chose p = 3500, hence = 1. Next, the parameter lu is given by 
u = tl>2 /wi, where lu\ and ui2 are the mean lifetimes of the attractiveness and of the active burglars, 
respectively. The values chosen in [18] are wi = (l/14)day _1 and 0)2 = 6 day -1 , hence a; = 84. The 
parameter tp is given by -0 = 9T /u>iu>2, where F, describing the source term, stands for the number 
of burglars becoming active per time and surface unit and 9 is the factor that boosts attractiveness 
due to a burglary. We take T — 0.002 burglars -day -1 - house separation -2 , which is the value taken 
in [22] for the analogous parameter (the value T = 0.5 used in [18] seems too high, apparently due 
to a lack of spatial normalization). For the domain that we consider, this means a total source 
term of 7 burglars ■ day -1 . Concerning 9, following [18], we take 9 — 1 burglars -1 • day -1 , which 
leads to ip = (14/3) x 10 -3 « 0.0047. As for the attractiveness diffusion parameter 77, which 
models neighboring effects, smaller values mean that repeat victimization is more and more likely 
than near-repeat victimization. The values taken in [22] range from 0.01 to 0.2, whereas [18] takes 
77 = 0.001. 

Last, let us consider the initial data. Concerning the (adimensional) attractiveness A, recall 
that the constant A is assumed to represent its static component, while A — A represents the 
dynamic component associated with the boost hypothesis. However it is not made completely 
clear in [22J [18] how A should be evaluated or measured. Here we shall assume that 

A < A {x) < 1, (10) 

which implies that A m i n = A and A max = 1. Consequently (cf. Lemma 12.11 below), the attrac- 
tiveness remains for all times larger or equal to A, and at most one. Note that in the numerical 
simulations from |18| . the initial data are assumed to be close to the homogeneous steady-state: 



/"ip - I + J (ip - I) 2 + AipA \ 

(^*)=( — ^ 4 
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which in particular satisfies (|10p. As for the initial burglar density No, we assume for simplicity 
that its average value N satisfies 

j^:=|nr 1 ||jv ||i<i = JV*. 

In dimensional units (cf. [TB]), in view of the above choice of 0J2 and T, this corresponds to an 
initial total population less than ps 1.2 burglars (although this may seem small, recall that we have 
at the same time a total source term of 7 burglars • day" 1 ). 

Now, for the above-chosen numerical values of the parameters, the global existence and bound- 
edness condition in Theorems 11.11 and 11.31 becomes 

A- 2 (l-A) < 7 7,, where 7 = * (|fi|^)-i = « 10.31. (11) 

This is equivalent to 

i>i_:= 2 



1 + y/l + 4777 

The values of A_ as a function of 77 are given in the following table: 



V 


A_ 


0.01 


0.91 


0.05 


0.73 


0.1 


0.61 


0.2 


0.49 



Therefore, our global existence conditions are compatible with magnitudes of attractivity which 
are up to about twice those of their static component in the range of rj used by [22] . We note that 
it is suggested in 18 that a ratio of order 10 (instead of 2) might be desirable, even for smaller 
values of r\ that the ones we consider here. However, this is still beyond the range in which we 
can rigorously prove global existence (and actually, global existence for such systems need not be 
taken for granted - see next subsection). 

1.3. Related results and comments. As noted in [22], problems ([I]) and (j3j belong to the 
family of chemotaxis-type systems. A more general class of such systems takes the form 

r,AA-A + f(A,N), 

(12) 

AiV - V(NVh(A)) -ujN + g(A, N). 

The typical feature of chemotaxis-type systems is the presence of the antidiffusion term, here 
— \7(N\7h(A)), by which the population, of density N, tends to move towards higher concentrations 
of A. The best-known among this class of problems is the Keller-Segel model [14], corresponding 
to / = N, h — A and g — 0. There is a very large literature on such systems (see [TTJ [T^], for 
recent surveys). 

System (|12p with a logarithmic sensivity function h = x l°g(A) was considered in [U [23] for 
f = N and g = lu = 0. Without any size restriction on the initial data, existence of a global classical 
solution was proved in [4] under the restriction ^ < 1, whereas existence of a (possibly singular) 
global weak solution was proved in [23j for any \ > 0. System (|12j) with a logarithmic sensivity 
function was also considered in [2] (see also |17j). in conjunction with a quadratic absorption 
term, namely g = —cN 2 , and / = N. In this case, global classical existence was proved for any 
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X > without size restriction on the initial data, but the proof made crucial use of the quadratic 
absorption term. 

Here, with only a linear absorption — ljN and the specific / under consideration, a size condition 
on the initial data is required to prove existence of global classical solution, but we do not need 
absolute restrictions on the parameter x- Although we do not know what happens for larger 
initial data, let us recall that in the case of the Keller-Segel model in two space dimensions, global 
existence is true only for small initial mass and that blow-up may occur for large mass. However 
proving blow-up for 2-dimensional parabolic chemotaxis models is a very challenging task. The 
only known blow-up result for the 2d Keller-Segel system is that of Herrero and Velazquez [TO] , 
and its proof, based on matched asymptotics, is highly technical. In simplified parabolic-elliptic 
chemotaxis models, where the term dA/dt in (fT2|) is replaced by 0, blow-up proofs are easier and 
there are more results available. In particular, when h = xlogA, f = N,g = uj = and the 
space dimension n is at least 3, finite time blow-up is known |15] to occur for x > 2n/(n — 2). For 
blow-up results concerning the parabolic-elliptic Keller-Segel system, see e.g. the survey [BJ. 

Our global existence proofs, as is customary in chemotaxis-type problems, rely on a priori 
estimates obtained by energy and entropy arguments. However, in the case of Theorem ll.il some 
care is needed to properly take into account the size condition on the initial data. In the case of 
the Keller-Segel system, an exact entropy functional was found and used to prove global existence 
for suitably small mass in [TBI [91 H]. Although system © is not known to possess an exact entropy 
functional, assuming condition ^ on the data of the problem, we can nevertheless construct an 
approximate entropy functional (with the typical N log N growth - cf. Lemma 13. 1[) , which is the 
key to our a priori estimates. We also stress that, in spite of uniform lower and upper bounds for 
A which easily follow from the maximum principle, the proof of the estimates of N is far from 
being immediate. As for Theorem II .21 under the assumptions % < 1, .A < 1 and nmxAo < 1, we 



can use a modified entropy functional from [5]. 

Remark 1.1. By straightforward modifications of the method (see Remark |4. II for more details), 
one can show that a result similar to Theorem 11.11 remains true for the initial-boundary value prob- 
lem associated with the more general system (|12[) where /, g, h are (sufficiently smooth) functions 
satisfying the following conditions 

• (Hi) g(A,0) > 0, VA > 0. 

• ("^2 ) There exist two positive constants A m ; n < A maj[ such that 



+ /(A min ,A0>0, 

n max 



ViV > 0, 
V7V > 0. 



\g(A,N)\ < gi(A)N 1 ^ s + g 2 (A), 
where 5 £ (0, 1) and g\,gi are bounded in L4 m ; n , -A max ]. 



• («4) 




sup 

1 „,A max ] 



d l h{A) 
dA l 



< oo, i 



1,2. 
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The outline of the rest of the article is as follows. Subsection 2.1 is devoted to local existence 
and uniqueness, whereas Subsection 2.2 contrains basic estimates, energy identities and functional 
inequalities which will needed for proving global existence. Sections 3-5 are then devoted to the 
proofs of Theorems 1.1, 1.3 and 1.2, respectively. Finally, our main conclusions are summarized in 
Section 6. 

2. Preliminaries 

2.1. Local existence and uniqueness. We have the following local existence and uniqueness 
result. 

Proposition 2.1. Assume fl|), and fix any (3 £ (0, |). Then problem admits a unique, 

classical, maximal in time solution (A,N) with A > 0. Moreover, if its maximal existence time T* 
is finite, then lim t ^ T » ||A(i)||#i+/3 + ||Ar(f)|| L 2 = oo. 

Note that the solution given by Proposition l2.1l is of course independent of the choice of (3 £ (0, |) 
(due to local uniqueness and the fact that ff 1+ ^(f2) c iJ 1+7 (fi) for (3 > 7). The procedure for 
proving Proposition ^. H is rather standard (see, e.g., [17j[2j[T3] for more details on similar problems). 
Although the function i9 is singular at A = 0, this actually causes no difficulty, because a positive 
lower bound for A can be deduced from the maximum principle (first working with a smooth 
replacement of We sketch the proof for completeness. 

We first recall some properties of theory of abstract evolution equation. Let W, V be two 
separable Hilbert spaces with dense and compact embedding V C %. Identifying H with its dual 
H' and denoting the dual space of V by V. 

Consider the initial value problem 

f +AU = F(U), 0<t<oo, (13) 
U(0) = U . 

Here, A is bounded linear operator from V to V which is defined by a symmetric bilinear form 
a(., .) on V satisfying 

|o(I7,V)|<M||l/||v||V||v s U,V€V, 

a{U,U) > c\\U\\l, UeV, (14) 

with some positive constants c, M. F : V — > V' is a continuous mapping such that, for each 6 > 0, 
there exist nondecreasing functions 4>g , ipg such that 

\\F{U)\\v> <0\\U\\ v + M\\U\\h), (15) 
\\F(U) - F(Y)\\v> <e\\U- V\\ v + (\\U\\v + \\V\\ V + 1)M\\U\\h + \\V\\h)\\U - V\\ n , (16) 
for U, V € V. Then by standard argument (see [20]), the following holds: 

Proposition 2.2. Assume |i^|)-( T76]) and let M > and Uo £ H with \\Uo\\n < M. Then there 
exists T = T(M) > and a unique local-in-time solution U to \13\) such that 

U € C([0,T];H) n i/ 1 ([0, T]\V') n L 2 ([0,T];V). (17) 

To deduce Proposition 12.11 from Proposition 12.21 we fix (3 in (0, |) and we set 

V := H 2 N +f3 (Q) x H\n), H := H 1+(s (n) x L 2 (fi). 

Then 

v' = H /3 (f7)x(H 1 (r2))'. 



8 



R. MANASEVICH, Q. H. PHAN, AND P. SOUPLET 



The linear operator A is defined by 

A = 



Ti 
T 2 



where T± = — A + 1 is regarded as an operator from Hj^~^(Q) to H^(fl) and T 2 = —A + ui is the 
Laplace operator equipped with the Neumann boundary condition in H l (VL)' . Then the bilinear 
form a(., .) is defined by 

a(U, V) = (T[ 1+P}/2 A 1 ,T^ +0)/2 A 2 ) + [ (V7V lV A 2 + uN^dx, 
v ' L2 Jn 

where U = (Ai,Ni) and V = (A 2 ,N 2 ). 

For any fixed a £ (0,A m - m ), we pick a function d a £ C°°(R) such that i? a (s) = xlogs feti- 
sh A mln . Define ^(C/) by 

F a (U)=( <PNA(1-A)+A \ 

and denote by (P a ) the modified problem (jHJ. 

As a consequence of Proposition 12. 2\ problem (P a ) admits a unique, maximal in time solution 
(A, N) , defined in an interval [0,T*). Moreover, ifT* < oo then lim t ^ T « \\A(t)\\ H i+f, +\\N(t)\\ L 2 = 
oo. Furthermore, it is classical for t > 0. This follows from a standard bootstrap argument based 
on parabolic regularity and imbedding theorems. 

Observing that A m ; n is a subsolution of the equation for A, we deduce that A > A min in (0, T*), 
so that (A, N) actually solves the original problem ([3]). 

As for local uniqueness for the notion of solution to ([3]) implies that A is uniformly positive 
onSJx [0, T] for any < T < T* (since A £ C([0, T*); H 1+0 (fl)). Since a can be arbitrary close to 
in problem (P a ), the uniqueness for the original problem follows from the uniqueness for (P a ). 

2.2. Basic estimates, energy identities and functional inequalities. Our first lemma pro- 
vides the primary L 1 control of N and uniform lower and upper bounds for A. The latter can be 
easily obtained from the maximum principle, owing to the special form of the nonlinear term in 
the equation for A. 

Lemma 2.1 (A priori estimates for A and N). For all t £ (0,T*), we have 

A(x,t) > min{l, A, mi A (x)} := A min , (18) 

A(x,t) < max{l,i,supA)(x)} := A max , (19) 

N(x,t) > 1 - e-^ > 0, (20) 

11^)11! = e-^||ATo||i + |0|(1 - *), (21) 



hence in particular 



|JV(f)||i < max{||JV ||i, \Sl\} := N 1>max . (22) 



Proof. To check (fT5|) and (fl9"]l. it suffices to observe that A min ,A max are, respectively, sub-/super- 
solution of the equation for A, with Neumann boundary conditions, and to apply the maximum 
principle. The proof of (|20|) is similar, in view of the regularity of the solution. As for (|21| . it 
follows by integrating the equation for N in space. □ 



We next state the basic energy identities. 
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Lemma 2.2 (Energy identities). For all t £ (0,T*) 7 we have 
1 d 
2dt 



li\\A(t)\\l + \\Mtm+r)\\VA(t)\\*=i> I NA 2 (1-A)dx + I AAdx, (2Hi 



~^l|VA(t)||l + ||VA(t)||l + v\\AA(t)\\l = -il>J NA(1 - A)AAdx, (24) 

/ (jVlogvV-iV+l)dx + w / (iVlogiV-iV + l)cfa; + / dx - [ VN.V$(A)dx 

dt Jo Jo Jo N J n 

= uj [ (log N — N + \)dx < 0, (25) 



^T-H iV (*)lli +ti; ll iV '(*)lli + ll ViV (*)lli = / NVN.V$(A)dx + oj [ Ndx. (26) 
2 ere Jo Jo 

Proof. Formulae (|23l) - (j2"6")l follow by integration by parts after multipying, respectively, the first 
equation in ^ by A and — AA, and the second equation in ([3]) by log N — 1 and N. Note that 
these manipulations are licit in view of (|20p. of the classical regularity of N and of the higher 
parabolic regularity applied to the equation for A. □ 

Lemma 2.3 (Poincare's inequality). For any uniformly positive function u G _ff 1 (J7), there holds 

Hl!</i 2 Hi / Mlds + ini^Hi?, (27) 

Jo u 

where fi is the best constant of the following Poincare-Sobolev inequality: 

\\u- |Oj— 1 / udx\\ 2 < Ml|Vu||i, m6W m (0). (28) 
Jo 

Proof. Denote u = L udi. Using the orthogonality of u — u and u in L 2 and (|28l) . we have 



lull 2 = \\u-u\\ 2 + \tt\u 2 



< /i 2 ||Vu|| 2 + \n\u 2 . (29) 
Since ||Vu|| 2 < ||u||i J n dx by the Cauchy-Schwarz inequality, we deduce (|27p. □ 

A key ingredient in the proof of our main result is the following interpolation estimate. We 
recall that H 2 (il) C C(f2), since we are in two space dimensions; however Lemma [2~4l remains true 
in any dimension if we assume u € H 2 ( fl) n (7(17). 

Lemma 2.4. Assume that u G H 2 (fl) satisfies du/dv = on dfl (in the sense of traces). Then 



[ l^u^dx < K osc 2 (w) / |Au| 2 da;, 
Jo Jo 



(30) 



where K = #(fi) > 0. 

Proof. Step 1. We first give a homogeneous version of the standard elliptic L 2 -estimate, namely: 

1 /2 

||D 2 u|| 2 := (^K^-lli) < C(fi)||Au|| a , (31) 

ij 

for any u € H 2 (Q) with Neumann boundary condition. 

To verify (|3"Tj) . we start from ||D 2 u|| 2 < (7(0)11 — Au + u\\ 2 , which is well known (see e.g. [7[ 
Chapter 9]), hence 

\\D 2 uh < C(n)(||A«|| a + |H| 2 ). (32) 
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Let (efe)fe>o be a Hilbert basis of L 2 made of eigenfunctions of —A with domain H 2 equipped 
with Neumann conditions. Denote the eigenvalues by Ao = < Ai < A2 < ... and observe that 
eo = C = Const. Set / = —Am, decompose / = ^2 k>0 c^efe, and note that 

c = (/, e ) = C I f dx = —C I d v uda = 0. 
Jn Jon 

Let v := X)fc>i 1<2 k£k- Then —Av = f with d u v = and v satisfies 

\\v\\ 2 2 =J2K 2 \ck\ 2 <\i 2 \\Au\\l (33) 

k>l 

Since the difference z = v — u satisfies Az = 0, d v z = 0, it follows that f n |Vz| 2 dx = 0, so that 

u = v + Const. 
Combining this with (|32|) and (|33[) . we obtain 

\\D 2 u\\ 2 = \\D 2 v\\ 2 <C(n)(\\Av\\ 2 + \\v\\ 2 ), 

and 

\\D 2 u\\ 2 <C(n)(l + X^)\\Au\\ 2 . 

Hence (HQ). _ 

Step 2. By density, it suffices to prove the Lemma for u S C 2 (f2). We may assume without 
loss of generality that minji u — and osc(it) = HuHoo. For a matrix M = (niij), we denote 
\M\i = maxij^j |mij|}. Observing that 

V ■ (|V?i| 2 Vu) = \Vu\ 2 Au + 2 t Vu(D 2 u)Vu = t Vu(2D 2 u + (Au)l)Vu, 
where / is the identity matrix, we have 

|V • (\Wu\ 2 Wu)\ < \2D 2 u + (A^ll^Wul 2 . 
Using the divergence theorem, it follows that 

/ \Vu\ 4 dx = / Vm • \Vu\ 2 Vudx = - / uV \\Vu\ 2 Vu)dx 
Jn Jn Jn 

< / u\2D 2 u + {Au)l\Ay U \ 2 dx <\ / |V?i| 4 dx+i / u 2 \2D 2 u + (Au)l\\dx. 
Jn 2 j n 2 j n 

Consequently, 

/ \Vu\ 4 dx < / \2D 2 u+ {Au)l\ 2 dx. (34) 

Jn Jn 

Since J n \2D 2 u+ (Au)l\^dx < C\\D 2 u\\l, it follows from ((3TJ) that 

/ |2D 2 u+ {Au)l\\dx < K{n)\\Au\\ 2 2 (35) 
Jn 

which, along with (|3"4"|) . implies the conclusion. □ 

We shall also need the following classical smoothing properties of the Neumann heat semigroup. 
For d, A > 0, we define the operator A = Ad.\ = —dA + A on L 2 (Q), with domain D(A) = {v £ 
H 2 (Q); dAjdv — on dtt (in the sense of traces)}. We denote by T(t) — Td,\{t) the semigroup 
generated by A. It is well known that, for any 0<m<2,l<p<<7<oo and any <j) G 
L 2 (n)nL p (fl), we have 

\\T{t)<t>\\w™.«{n) <Ce- xt (l + t-^-^)U\\ LPm , t>0, (36) 
with C = C(f2, m, p, q, 77) > 0. 



CRIME MODELING 



11 



3. Global existence and boundedness: Proof of Theorem [LT 



The following lemma is the key to the proof of our main result. It provides an approximate 
entropy functional <p, which is available whenever condition ([7]) is satisfied. This function enables 
us to get an 7J 1 -bound for A and a time-averaged 7J 2 -bound. 

Lemma 3.1 (Approximate entropy functional). Assume ^ and let 

ci ■■= \ (i - ^xVVVJv?,™A^(^ - A ^ 2 ) > °- ( 37 ) 

(i) Then the function 

f 1 2ib 2 

<f>(t):=a {NlogN-N + l)dx + -\\VA(t)\\l with a = ^-A^ 2 N hmax , 
Jn 2 r] 

satisfies the differential inequality 

(/}' + + Cl \\AA\\j < c 2 , 0<i<T*, (38) 
where uj = min(cj,2) and c 2 = c 2 (K, x, A min , A max , ip, JVi, max , r), |0|,/j) > 0. 

sup \\VA(t)\\l < c 3 , (39) 
te(o,T*) 

with C3 := 2 max(0(O), C2<2>~ 1 ), and 

/ ||AA(r)|||dr < c 4 (l + t-s), < s < t < T* , (40) 

with C4 = cj" 1 max(c2, C3/2). 
Proof. Set 

01 (t) := ||VA(t)||l, <j> 2 {t) := [ (N\ogN -N+l)dx>0 (41) 



(due to slogs — s + 1 > 0, s > 0). In this proof, C will denote a generic positive constant depending 
only on K, x, A min , A max , ip, Ni.max, »?, and on £i,e 2 ,e 3 below. 

On the one hand, it follows from (|24l) that, for any e\ > 0, 

/,2 



Since 



h[ + 0x + ^AA^ < ei ||AA||l + N 2 A 2 (l A) 2 dx. (42) 

||iV||| < ^ 2 iVi, m ax / ^r-dx + C (43) 



due to Lemma |2~31 we deduce from (|4"2"j) that 



where 



Vi +<t>i<ax ( ^-dx + a 2 \\AA\\ 2 2 + C (44) 
z Jn ^ 

01 = - — ^ ax /i 2 iVi imax > 0, a 2 = £1 - ?/. (45) 
4ei 



On the other hand, it follows from (|25j) that, for any e 2 , £3 > 



/> 2 + UJ(j) 2 



N Jn Jn N 4e 2 Jq 

s " /„ h^" 1 + i( E » l|Ar|li + i /„ |v " ( - 4)|4dI ) ' 
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Using 



/ |Vi?(A)| 4 ^ < X x 4 A m f n (A max - A Blin ) 2 \\AA\\ 
due to Lemma l2~4l and (143)) , we deduce that 



<>, -^.>> ' <':; I ^^dx + a 4 \\AA\\ 2 2 + C, (46) 



£2 



where 



a 3 =s 2 + 1, C4 = A Aj n (A max - ^ min ) 2 > 0. (47) 

Now setting uj = min(cj, 2) and combining ([33)1 and (|4"6")) . we see that = \<pi + o-(j) 2 satisfies 

r IViVl 2 

<f>' +uHf><(a 1 + aa 3 ) ^—^dx+(a 2 +aa 4 )\\AA\\ 2 2 + C. (48) 

Assume 03 < and choose a = —01/03 > 0. Then we have a 2 + oa 4 < provided 0104 < 0203, 
that is 

^X^'fn^fnax^max " ^mi„) W^lmax < lfexfa - £i)(4(l - £ 2 )£ 2 " ^M^l.max) £ 3 • (49) 

The best condition, maximizing the RHS in ([4TJ]). is obtained by choosing e\ — n/2, e 2 = 1/2 and 
then e 3 = (2/i 2 A f l max ) _1 , which in turn implies a 3 = —1/4 < 0. Inequality (|49fl is then equivalent 
to 

^^jn^Lf^x - A min )W 2 ^l,max < 4ry 2 (l - e 3 M 2 ^l,max)£3 = ^(/X^l.max)" 1 , 

which is true, due to (|37p . Then we have 02 = —i}/2, 

_ ^ 2 Ai 2 at _ -^X 4 2 it ,.-4 / A a \2 _ ■^V-' 2 A 4 2 AT 

We conclude from flU} that f[38]) holds with ci = — 02 — 0-04, which yields the value given in (|37|) . 
(ii) Multiplying (l38l) with e"* and integrating between and t, we obtain 

cj)(t) < max(^(0), C2W" 1 ), < t < T*. (50) 

This guarantees ([39]) , in view of (|4T|) . Inequality (|40)) then follows after integrating (f38|) over (s, i), 
taking (|50|) into account. □ 



Building on estimate (|40|) from the previous lemma, we shall now derive uniform estimates for 
||7V(i)||2 and ||A(t)||#WQ\, which in turn will guarantee the global existence of the solution. 

Lemma 3.2. Assume that |^0[ ) is satisfied for some C4 > 0. Then we have 

sup \\N(t)\\ 2 < 00 (51) 

te(o,T*) 

and, /or eoc/i to € (0, 2) and r g (0, T*), 

sup \\A(t)\\ B m(o)<°°- (52) 

te(r,T*) 



CRIME MODELING 



13 



Proof. Let us first establish (|5"Tj). It follows from (|26|) that 

^||iV"(t)||l + 2a;||iV||| + 2||ViV(t)||| < 2uN 1]max + 2 / XXX.Y,H Ahi.nli 



<2uN Umax + / VN 2 .V$(A)dx 



i a 
vN 2 . 

>a 

<2wAT 1;max - / 7V 2 .Atf(A)dx 
in 

<2wJVi. max + ||JV||l||Atf(i4)|| 2 . 
By the Sobolev embedding of W 1,1 (f2) into L 2 (f2) with constant Cg, we have, for any E\ > 0, 

11^112 - \W% < c s (||viv 2 ||i + ||tv 2 ||i) = c s (2||jwiv|| x + ||jyr||l) 

< C s (2||iV|| 2 ||ViV|| 2 + ||iV|| 2 ) < e!||ViV|| 2 + (C s + C 2 / £l )||7V|| 2 . 
Choosing ei — ||A$(A)|| 2 , we get 

|||7V(t)|| 2 + 2u;||iV|| 2 + || W(i)|| 2 < 2 w JV 1;max + (C g \\A0{A)\\ 2 + C 2 ||Atf(A)|| 2 ) ||7V|| 2 . (53) 

Pick now e > 0, and denote by C(e) a generic positive constant depending on the solution and 
on £, but independent of t € (0, T*). By the Gagliardo-Nirenberg inequality 

||iV|| 2 <C G (||V7V||^ 2 ||iV||l /2 + HJVIU), 

we have 

||7V|| 2 < 2C 2 (\\VN\\ 2 \\N\U + ||7V|| 2 ) < £ ||ViV|| 2 + C( £ )7V 2 max . 
Therefore, using (]22f) . we have 



HViV|| 2 >i||7V|| 2 -C( £ ). 



Combining this with (|53|) . we obtain 
Now setting 



\\N(t)\\i + -\\N\\i< C(£) + (1 + 2C 2 || A^)|| 2 )||iV|| 2 . 



p(t) = ||iV(t)|| 2 , P {t) = jT*(i - 1 - 2C 2 ||Atf(A)|| 2 )d s , 
we are thus left with the differential inequality 

vt{t)+ff{t)<p{t)<C{e). 
On the other hand, we have Ai?(j4) = x{A^AA - A~ 2 \X/ 'A\ 2 ) , hence 

||A 1 ?(A)|| 2 <2 X 2 (A- 2 ||AA|| 2 + A- 4 ||VA||4). 
It then follows from ([ID]). (TT51) and Lemma EH that 

f ||Atf(A)|| 2 <c 5 (1+*-s), 0<s<t<T*. 

By choosing e = (2 + 2CgC5) _1 and letting cq = 2C|c5, we have 

p{t) - p{s) > t - s - c 6 , 0<s<<<T*. 

By integration, we get 

tp(t) < <p{0)e~ p ^ AC I e p(s) - p{t) ds < ¥ j(0)e C6 "* + C / e C6+s ^ds < e c %(0) + C) 
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and (|51~j) follows. 

To prove (|52|). we rewrite the first equation via the variation-of-constants formula 

A(t) = T., hl {t)A Q + f T ritl (t-s)(i>NA(l-A)+A)(s)ds 
Jo 

(where T„ ; i is defined at the end of Section 2). By (|51l) . we have 

M:= sup \\(ipNA(l -A) +A){t)\\ 2 < oo. 

t£(0,T») 

Fix to £ (0, 2). It follows from ([36j that, for any < t < t < T* , 

\\A(t)\\ Hm{n) < Ce-\l + <- m/2 )|| A || 2 + C / e-(*- 5 )(l + (t - .s)- m / 2 )|| (^AM(1 - A) + A)(s)\\ 2 ds 

Jo 

poo 

< C(l + r- m / 2 )||A || 2 + CM / e" s (l + S - m / 2 )ds =: C(r), 

Jo 

which proves ([52")) . □ 
Proof of Theorem 

Step 1. Global existence. In view of the local theory stated in Section 2, this is a direct 
consequence of estimates (fSTj) and (|52"|) in Lemma 13.21 

Step 2. Boundedness. Estimate ([5]) was already obtained in Lemma 13.21 Starting from the 
global estimates obtained in Lemma [3~2l we shall use a standard boostrap argument (see e.g., |16j ) 
to prove ©. 

As a consequence of (|52|) and Sobolev imbeddings, we have, for each t > 0, 

sup ||VA(t)||p < oo, \<p<oo. (54) 

te(r,oo) 

Using (fTHl) and ((5TJ) , it follows that 

M q := sup IKiVA^VAX*)!!, < oo, l<g<2. (55) 

te(r,oo) 

Set T(t) = T\ yUJ (t). From the second equation, using T(£)l = e~ ut , we have 

N(t)=T(t)N + [ T{t- s)(-V(NVtf(A)))ds+ [ T(t-s)u)ds 
Jo Jo 

= T(t)N + 1 - e-"* + [ VT(t-s)-NA- 1 VAds. 
Jo 

After a time-shift, combining this with (|3"6")l . (|5ip and (|55j) . we obtain, for any e [1, 2), m € [0, 1) 

and t > t/2, 

||^(f + f)lk~.«(n) < Ce-\l + t- m/2 )||iV(|)||« + C 

+ <?/ e-( t - s )(l + (t-s)- I2 ^)||(ArA- 1 VA)(| + s)|| (? ds (56) 
Jo 

< C(l + r-™/ 2 ) + CMg / e- s (l + s-(" i+1 >/ 2 )ds =: <5(r). 
Jo 

Using Sobolev's imbedding again, taking q close to 2~ and to close to 1~, we deduce 

sup ||AT(£)|| p < oo, 1 < p < oo. 
ie(r,oo) 
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Due to (154]) . this implies that (|55|) is true for any q £ [l,oo). The argument in (|56|) then yields 

sup ||-AT(t)||vf m .«(n) <°°> < m < 1, 1 < g < oo 

te(T,oo) 

and ((9]) follows by a further application of Sobolev imbeddings. □ 
4. The case of a square domain: Proof of Theorem 11.31 

Let us consider the square f2 = (0,L) 2 . Theorem 11.31 is a consequence of Theorem 11.11 and of 
the following two lemmas, where we estimate the constants K and fj, in inequalities (|30[) and (|28p 

Lemma 4.1. for f2 = (0, L) 2 , the constant K in estimate UJfy) satisfies K < 12. 

Proof. Recall that the constant K = K(il) is determined from inequality (l3"5j) in the proof of 
Lemma 12.41 Repeatedly using the Cauchy-Schwarz inequality, we first compute 

\2D 2 u + (Au)l\ 1 = max{|3 

= 2\u xy \ + max{|3 

H-xx ^yy 1 7 | ^xx <^^yy \ | f 

< 2| U:C!/ | + (io(t& + 1/2 

and then 

\2D 2 u+(Au)l\l < 12(u 2 xx + u 2 yy + 2u 2 xy ). 

Consequently, 

/ \2D 2 u+(Au)l\ldx<12(\\u xx \\l + \\u yy \\ 2 2 + 2\\u xy \\ 2 2 ). (57) 
Ja 



We next claim that 



KxB + IK3/II2 + 2\\u xy \\i = \\Au\\i (58) 



Indeed, using the Fourier expansion u = ^2- k>Q ajk cos cos -p^, we have 



l xx 1 1 2 "i - II '"yy II 2 



^E^o + ^ E fc4a o fe + ^ E 0' 4 + fc4 K 2 ^ 

j>l k>l j,k>l 



^4 

will = E j 2k2a %> 

j,fc>l 



^EA 2 » + ^ E fc4 <4 + ^ E O' 2 + fc2 ) 2 <4> 

j>i fc>i i,fc>i 

which implies (|58p. 

The conclusion follows by combining (f3"5"j) . (|57| and (|58|l. □ 

Lemma 4.2. For fl = (0, L) 2 ; the constant [i in estimate (28\) satisfies /i < 

Proof. Let u be a smooth function with zero average. There exist mean values C,{x) 1 ^{y) in [0,L] 
such that 

1 f L 1 f L 

i(s, y)ds = u(£(y),y), — \ u(x, t)dt = u(x, ((x)). 



It follows that 

u(x,y) = Y I u(s,y)ds+ [ u x (s,y)ds 



LJ 
1 

I 



£(») 

l rv 
u(x,t)dt+ / u y {x,t)dt. 

k(x) 
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Consequently, 

u 2 (x,y)dxdy =J2^J u(x,y)dxdy 

+ I (— J u(s,y)ds I u y (x,t)dt\dxdy 
Jn \L Jo Jcfx) J 



l_ ' 

L 2 



(— I u(x,t)dt j u x {s,y)ds\dxdy 

u x (s,y)ds / u y (x,t)dt \dxdy. 
Q \J£(y) JC(x) " ) 



Since u has zero average, it follows that 



1 



IM| 2 < IMI 1 (|| 1t a ||l + ||%||l) + K||lK||l. 

Moreover, by [1] Theorem 3.2], we have ||u||i < (i/v2)|| Vit||i, hence 

||u||i < ^||Vu||i(||ux||i + K||i) + \iWuxWx + Kill) 2 . 

Using the elementary inequality ||ux||i + < v2||Vii||i, we deduce that 

||u|| a <y|||Vu||i. 

Since C x {9) is dense in W 1 ' 1 ^), Lemma I4T21 follows by letting u = N- J n N(x,y)dxdy. □ 

Remark 4.1. Let us briefly justify the generalization in Remark 1.2. Conditions ("Hi), (H2), (H5) 
are natural. By condition (K.3) and N x ~ s < eN + C(e), one can show that ||JV(£)||i is uniformly 
bounded. And then, one has following estimates 

| 5 (4JV)logJV|<Ci(A)^ + C 2 (A) (a) 

\g(A,N)\N<eN 2 + C 3 (e,A) (b) 

Inequality (a) and boundedness of ||iV(i)||i imply J46| ). Inequality (b) implies \5S\) . By condition 
{T-La), one has 

\f(A,N)\ 2 <C 4 (^)7V 2 + C 5 (^) 

which implies inequality 

5. The case x < 1: Proof of Theorem 11.21 
For x < 1) a basic entropy estimate is given by the following. 
Lemma 5.1. Assume < \ < 1 and j4 max = 1- Then 

sup / N\\ogN\dx < 00. (59) 
o<t<T* Jn 

Proof. Following [5J Section 4], we define 

Y(t)= / N(\ogN - clogA)dx. 
Jn 
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Letting F = ipNA(l - A) + A - A and G = u(l- N), we have 

Y'(t)= [ N t (logN -c\ogA)dx+ / (JVf — cNA~ l A t ) dx 
Jn Jn 

= Q+ I (log N-c log A)Gdx + [ Gdx-c [ NA~ 1 F dx, (60) 
Jn Jn Jn 



where 



Q = [ \aN-W • (JVVi?(A))1 (log TV -clog A) dx~cn [ NA' x AAdx 
Jn L J j n 

N- 1 \WN\ 2 dx+ (x + c(l +77)) / A _1 ViV- VA-c(x + ry) / iV^" 2 |VA| 2 rfx 
a Jn Jn 



= - / N-'A- 2 
Jn 



A 2 \WN\ 2 + c(x + Ty)A^ 2 |VA| 2 - (x+ c(l + ?7))iV^ViV • VA 



We see that Q < provided the discriminant <5 := (\ + c(l + r?)) 2 — 4c(x + 77) is nonpositive, that 
is 

(1 + r]) 2 c 2 - 2c{ X + 277 - xv) + X 2 < 0. (61) 
There exists a constant c satisfying (I6ip if and only if 

\x + 2v-xv\ > 0- + v)x> 

which is equivalent to x < 1. 

On the other hand, we have G < uj and, since A < 1 due to (JT5J), F > -A. Using ([TS]). (j22l) 
and the inequality log s < s (s > 0) , we deduce from ([60)) that 



y'(fj) < -ujY + lj / (log7V-clogA)d:r + w|ft| + c / iVdx< -ojY + C. 
Jn Jn 

We infer that sup 0<t<T , Y(t) < 00, so that ([59]) follows from l[l9"]). ([22]) . □ 

With the bound ()5.1[) at hand, we can now use the following better Sobolev type inequality in 
place of Lemma [ 



Lemma 5.2. For any e > 0, there holds 

HI < e I u\\ogu\dx [ LYHLdar + c'(n > ||u|| 1 ,e) (62) 
Jn Jn u 

for any uniformly positive function u £ H 1 ^). 

Proof. Fix > 1. Since V[u — fc] + = X{u>k}^ u , we deduce from (|29f that 

/ [u - k]\dx < n 2 ( [ \Wu\dx] + ( f[u-k] + dx 

Jn \Ju>k ) \Jn 
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On the other hand, we have u 2 < [u — k]\ + 2ku (just consider separately the cases u < or > k). 
Therefore, by integration, we obtain 

iril.r ; it- I / Vu\dx) +\n\- 1 \\u\\ 2 1 +2k\\u\\ 1 



[ 



O \Ju>k 
2 



, ;/ 2 i \jM_ dx . f U dx + \n\- 1 \\u\\ 2 1 +2k\\u\\ 1 

u>k u Ju>k 



(63) 



r IV7/I 2 / u \ lo & u \ dx 

< — l -dx- Jn - + \n\- 1 \\u\\ 2 + 2k\\u\\ 1 

Jn u log k 

and the Lemma follows by taking k = exp(/i 2 /e). □ 
Proof of Theorem \1.2\ We claim that 

sup ||VA(t)||| < oo (64) 
te(o,T*) 

and 

rt 

2 



I 



\\AA{T)\\i.dr <c(l+t-s), 0<s<t<T*. (65) 
Fix e > 0. By Lemmas 15.11 and 15.21 we have 



||iV|||< £ / W-dx + C(e), 0<t<T*. (66) 

Using estimate (l66l) instead of (|27|i in the proof of Lemma f3.1f i). we obtain the differential in- 
equality (|4"5)l where the constants dj are defined by formulae (l4"5"j) and (I4T1) with ^t 2 A r l max replaced 
by e. Like before, we assume 0,2, a 3 < and choose a — —ai/a 3 > 0. Then we have et 2 + 1704 < 
provided 0104 < 0203, which is now equivalent to 

^X 4 ^mf„^ma X (^max " Ani„)Ve < ~ 6l)(4(l - £2^2 - e 3 e)e 3 . 

Choosing £1 = 77/2, e 2 = 1/2, e 3 = l/(2e), the condition becomes Kx 4 A^ 4 m A^(A nll!lx - 
A m in) 2 4' 2£ < V 2£ l which is verified for e suitably small. Estimates (|64l) - (|6"5")) then follow as 
in the proof of Lemma |3 . lf ii) . 

The rest of the proof of the theorem then relies on Lemma 13.21 similarly as in the proof of 
Theorem O □ 

6. Conclusion 

In this article, we have considered a nonlinear, strongly coupled, parabolic system arising in the 
modeling of burglary in residential areas. The system involves two spatio-temporal unknowns: the 
attractivity value of the property and the criminal density. The system is of chemotaxis-type and 
involves a logarithmic sensivity function and specific interaction and relaxation terms. This model 
has appeared in [18] , as a modification of the model of Short et al. [22] . In [18] some results about 
linearized stability/instability of the homogeneous steady-states and some numerical simulations 
suggesting the existence of hotspots were given. However, the (local and) global existence of 
solutions was left open. 

In this article, under suitable assumptions on the data of the problem, we have given a rigorous 
proof of the existence of a global and bounded, classical solution, thereby solving the problem 
left open in [18) . In the range of anti-diffusion parameter \ relevant for the criminological model, 
our sufficient condition for global existence roughly says that, at the initial time, the product of 
the oscillation of the attractivity value and of the total criminal population should not be two 
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large. Our proofs are based on the construction of approximate entropies and on the use of various 
functional inequalities. 

We have also provided explicit numerical conditions for global existence when the domain is 
a square, including concrete cases involving values of the parameters which are expected to be 
physically relevant. In such cases, in the range of diffusion parameter r\ used by [22], our global 
existence conditions are compatible with magnitudes of attractivity which are up to about twice 
those of their static component. It is suggested in [18] that a ratio of order 10 (instead of 2) might 
be desirable, even for smaller values of r\ that the ones we consider here. However, this is still 
beyond the range in which we can rigorously prove global existence (and actually, global existence 
for such systems need not be taken for granted, as shown by the existing chemotaxis literature). 
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